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Next-to- leading order correction to the one-particle inclusive cross section in the 
framework of high energy factorization is calculated. Numerical results for 
midrapidity region are compared with predictions of conventional calculations 
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1 Introduction 



The recently completed calculation of the NLO correction to BFKL Pomeron [|IJ 
involves as its ingredients the formulae for the cross section of tree level two- 
particle production and one-loop virtual correction to one-particle production in 
Quasi Multi Regge Kinematics (QMRK). These results can be used in studying 
the properties of the inclusive two-particle production at leading order 0, |3|] and 
of the one-particle inclusive cross section in the next-to-leading order at high 
energies. 

The total cross section in the next-to-leading order(NLO BFKL) was cal- 
culated in |1| and is now under an extensive discussion j|, The total cross 
section has a complicated structure due to delicate cancellations of singularities 
in various virtual contributions and infrared divergencies arising in integrations 
over parameters of real correction. Cancellation between divergent parts of real 
and virtual contributions is also needed for computation of the next-to-leading 
order one-particle inclusive production cross section, but here the structure of 
this cancellation is simpler and more transparent than in the case of total cross 
section. 

The problem of one-particle inclusive production in high energy hadron colli- 
sion was investigated, in the leading order (LO), in a number of works ||. The 
purpose of this paper is to calculate one-particle inclusive cross section to the 
next-to-leading order at high energies proceeding as far as possible with analyti- 
cal calculations and then turning to numerical estimates. 

The outline of the paper is as follows. 

In the second section we briefly review the results on the particle production 
cross sections obtained within the high energy factorization scheme and compare 
them to the expressions obtained using collinear factorization. 

In the third section we describe a calculation leading to the explicit expression 
for the one-particle inclusive production cross section in the next-to-leading order. 

In Section 4 some numerical results on particle production in central rapidity 
region are presented and discussed. 

Section 5 contains a brief conclusion. 

2 Particle production at high energies 

From the theoretical viewpoint the hadron scattering at very high energies is 
special in the way the hard degrees of freedom (partons) are formed from colliding 
hadrons. When the ratio of hardness of the process k± (which is the transverse 
momentum of produced particle) to the invariant energy of colliding hadrons, 
\/~S, is not too small (up to 1CT 2 ) it is possible to describe the structure functions 
of hadrons by taking into account only processes that contribute logarithms of 
kj_/-^QCD at leading order, i.e. by resummation of the a™ \n n (k±/ Aqcd) terms. 
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Such structure function is given DGLAP evolution equation 0. 

Since the emergence of the combination a s \n(kx/ Aqcd) implies strong or- 
dering of emitted particles in their transverse momenta up until a hard collision 
block, the transverse momentum of detected particle k± is parametrically bigger 
than that of any parton involved in the process. Therefore it is possible to calcu- 
late the cross section of the hard process using the initial on-shell partons. This 
prescription is known as collinear factorization [0] and leads to the well-known 
result for the production rate: 

da f . „ , ,n,da, 



dk 2 dyi 



J dy 2 x 1 f a (x 1 ,k 2 )-^-x 2 f b {x2,k 2 ), (1) 



where f a (x,k 2 ) is a structure function for the parton of type a and x\p = 
k ± {e^ + e^yVS. 

At high energies for particles produced with k± <C another big logarithm, 
ln(l/x), is important. The resummation of such logarithmic contributions can 
become more important than of \n(k 2 /Aq CD ). The resummation of the leading 
energy logarithms for the structure function is described by BFKL equation ||. 
The domain of validity of the BFKL equation in describing the structure functions 
is at present not well understood. It is likely that for some kinematical region the 
correct approach is to resum the logarithms of both types, or at least interpolate 
between two types of resummation as done, e.g., in the CCFM|§ equation. 

The main point is that at high energies the transverse momenta of the incom- 
ing parton fluxes can no longer be neglected. To take them into account a new 
approach called k± or high energy factorization was proposed [|1(], ITf . Extensive 



description of the method and various applications can be found in fTT|| . Let us 
note that this method was de facto used earlier in [T2| . 



The method of high energy factorization is based on consideration of "partons" 
with nonzero transverse momentum that are, in contrast with the traditional 
collinear factorization case, virtual particles. In this case colliding hadrons are 
described by unintegrated structure function, 0, so that 



2 dxg(x,q 



0(^) = g 2 ^^, (2) 

where <fi/q 2 is proportional to the probability to find the incident parton (essen- 
tially a gluon, because gluons give leading contribution at high energies) with 
the longitudinal momentum component xp a (p a is a momentum of the incident 
particle) and transverse momentum component q±. 

Note that such an interpretation is somewhat oversimplified, because due to 
the quantum structure of QCD evolution some form-factors may arise changing 
the form of <ft. For DGLAP evolution it is a Sudakov form- factor []nj. However, 
when studying the semi-inclusive quantity like one-particle production cross sec- 
tion it is legitimate to use the unintegrated structure function in the simple form 
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ofEq.(0)Q. 

Scattering of the off-shell "partons" are described by generalized cross sections 
calculated in Quasi Multi Regge Kinematics (QMRK) fl4] , 15]. 



In QMRK one studies the 2 — > n + 2 scattering process with two outgoing 
particles having almost the same momenta as the incident ones and remaining n 
particles emitted into the central rapidity region separated by large rapidity gaps 
from incoming particles (the situation, where rapidity gaps between n particles 
are also large, corresponds to Multi Regge Kinematics, MRK). Large rapidity 
gaps allow to distinguish the quantities related to the incident particles from 
those describing the cross section of the hard process of particle production in 
the central rapidity region. 

Let us for example consider the cross section of the process gg — > ggg in the 
limit of high energy: 

dtJgg^ggg 4iV|ag t (ggl_L S {2) (gl_L + ?2_L ~ k±) ,„x 

d 2 k ± dy ~ tt 2 (N 2 -1) J q 2 ± k\ q 2 ± ' 1 j 

where qx t2 = p a ,b — p' a ,b- i n Eq. (§) the above-mentioned factorization of the 
cross section is clearly seen. Indeed, the first, second and third factors under 
the integral correspond to p a — > p' a ,qi splitting, (71,(72 ~~ > k "scattering" and 
Pt P21 12 splitting respectively. 

The factors related to the splitting of the incident particles should further be 
transformed to structure functions. This can be done in two steps. First, one 
assembles incident partons into wave packets describing by form factors || . The 
second step is taking into account additional radiation and virtual corrections 
summed to give the unintegrated structure functions <p(x, q±) with x fixed by the 
kinematics of the considered process. 

The cross sections of producing n = 1, 2 particles in the central region to the 
lowest perturbative order read: 

rf(J i f J2 j2 <p{xi,o, Qi±) dai <p(x 2 , , q 2 ±) 



J d 2 q 1± d 2 q 2 ± 



d 2 k±dy J ~ qf ± dk\ q 2± 

d&i 4:N c a s S^(q 1± + q 2± - fejj 
dk\ ~ N 2 — 1 k\ 

£1,0 = k ± e y /\/S, x 2 ,o = k±e~ y /vS] 
da 2 f 2 2 if(x 1 ,q 1± ) da 2 p(x 2 ,q 2± ) 



(4) 



J d 2 q 1 _ L d 2 q 21 _ 



d 2 k 1 ^d 2 k 2 ±dy 1 dy 2 J "~ " qf ± d 2 k 1 ^d 2 k 2 ±dAy q\ L 
da 2 2N 2 a 2 6^ (q 1± + q 2± - k 1± - k 2A _) 



d 2 k 1± d 2 k 2± dAy (N 2 - l)n 2 q 2 ± q 2 ± 



A, (5) 



2 I am grateful to Yu.L.Dokshitzer for pointing me this out 
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xi = k 1± e yi (l + k 2X e Ay )/VS, x 2 = k 1L e~ yi (l + k 2± e~ Ay )/^S. 

Ay = y 2 - yi- 

A part of the analytical expression for A can be found in flTJ, [Tj| for the 
subprocesses gg — > gg and in |15|, [17j for the gg — > qq ones. The explicit form of 



A was recently derived in [18| and is given in Appendix A. Note that the formula 



for gg — > qq cross section from |L7j coincides with the analogous formula in |TT 
in the limit of massless quarks. 

Eq.(f|) gives the rate of one-particle production in the leading order. It was 
studied in a number of publications [||]. Our aim is to calculate first correction 
to it. 

3 Real and virtual contributions to NLO one- 
particle production 

The one-particle production in the next-to-leading order includes two contribu- 
tions, real and virtual. The real contribution comes from the two-particle cross 
section Eq. (H) integrated over the phase space of one of the particles (considered 
unobservable) and the fixed four-momentum of the second particle: 

d °r _ f d ^ <p(xi,q 1± ) da 2 ip(x 2 ,q 2± ) 

j~f2i ,-7*2 7,. ,-7*2 7.1 s3 A n . ™ 2 * • 



d 2 k±±dyi J qf ± d 2 ki±d 2 k 2 ±dAy q% ± 

where 

rf$ = d 2 q 1± d 2 q 2± d 2 k 2± dAy (7) 

The factor of 2 in Eq. (§) reflects the identity of outgoing particles. 

The virtual contribution has the same form as in Eq.(^), but instead of &i we 
must use 

da v 4:N c a s 5 ( - 2) (q 1 _ L + q 2± - k ± ) 

dkj = Nf^i ¥ ± V(gi± ' q2±) (8) 

The virtual correction contains both ultraviolet and infrared divergencies. The 
ultraviolet one leads, through standard renormalization procedure, to the running 
coupling constant. The infrared divergence must cancel with the infrared and 
collinear divergencies in the real contribution. 

3.1 Cancellation of collinear and infrared divergencies 

Let us now outline at the formal level how this cancellation occurs. To deal with 
the collinear singularity we must introduce a jet defining algorithm which in two 
particle production case could be expressed through the function S(k, fci, k 2 ) {k\, 
k 2 and k are on-shell 4- vectors) so that Eq.(|6|) is replaced by 

A2 ? r A = I <Pk 1 ±d i k2±dyidy2-~ — , , S(k,ki,k 2 ). (9) 
d z k±dy J d z ki±d z k 2 ±dyidy 2 
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To provide the sought for cancellation between the real and virtual corrections 
S should be an infrared safe quantity, which means that the following property 
should hold (cf. 



3{k, Xh, (1 - X)h) = 6W(k x - h ± )5(y - Vl ), < A < 1 (10) 
In the following we choose S in the following form: 

S(k, ki, k 2 ) = 9(R>R )J2^ 2) (k ± -k i± )5(y- yi ) (11) 



i=l,2 



+ 6{R < /.'„) 'V 2 i/.-. - k 1± - k 2± )5 (y - - In 1 ' " h±< 



2 fci ± e-w + k 2 ±e~ 



i/2 



where R 2 = (<fii — 4> 2 ) 2 + (Vi — V^) 2 ■ Although the last line in Eq . (|TT|) may look ar- 
tificial, it has the natural meaning. If we claim that " combined" particle, formed 
by two particles indistinguishable under given resolution, has definite rapidity 
and that 4-momenta of two particles are added up to form the 4-momentum of 
"combined" particle, we immediately arrive at Eq. (|Tl~D . It is straightforward to 
check that Eq.flTTD satisfies Eq. fllCf) . 

Grouping together Eqs.(|),(^|) and ([□]) we find: 



d 2 k±dy J qfj_ d 2 k±d 2 k 2 ±dAy q\^ 

+ f M ?&p±L ^ 7K %^V <*.), (12) 

J q( ± d 2 h±d 2 k 2 ±dAy q 2± 



where k%± — k± — k 2 ± and Xx t2 = e ±2/ y S/S* with £ = k\ + k\ + 2kik 2 ch(Ay) (see 
Appendix A), and ki = \ki±\ which is a notation we shall use from now on. Note 
that when R — > x\ — > x± t o and x 2 — *> x 2j o- Let us now rewrite the second term 
in Eq.(H) as: 

,( p{xi,qi±) (p(x 2 ,q 21 _) <p(xi, , qix) (f(x2,o, Q2±)\ da 2 



Qi± q 2 ± qi± ) d 2 k i± d 2 k 2± dAy 

^(^1,0) Qi±) da 2 <p{x2,o,q2x) 



,-», ,,, ; % ■■> 6(R<R )- (13) 

5J_l d z ki±d z k 2 ±dAy q 2± 

The integration in the first line is free from divergencies, while in the second line 
integrals over d 2 k 2 ± and dAy do not involve structure functions and could be 
done analytically. Before doing this integration we make a replacement 8(R < 
Rq) = 1 — 9(R > Rq) in the last term and then substitute Eq.flilf) to Eq. dl2|) to 
yield: 

da r 
d 2 k±dy 
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/ d q 1± d q 2± - 2 - 2 / d k 2 ±dAy 



Qi± ?I± J d 2 ki±d 2 k 2 ±dAy 

M L «*<*A f> ^^ e(R ( kM)>R») (14) 

\ qf ± d 2 k ± d 2 k 2± dAy q^ ± 



<f(xi, ,qi±) da 2 <p(x2,o, ?2±) 



d 2 ki±d 2 k 2 j_dAy q% ± 



qf± ql± ql± ?i± 



f/ " 2 < /,>, 



d 2 ki±d 2 k 2 ±dAy 



where we indicate that R(k,k 2 ) = ((y — y 2 ) 2 + (0 — 02) 2 ) 1//2 and R(ki,k 2 ) = 
{{yi ~ IJ2) 2 + (0i — 02) 2 ) 1//2 are different in the third and forth lines of Eq.flTil). 
The third line in Eq.flHj) has an infrared singularity when k 2 — > and the forth 
one has singularities when k 2 — > or k\ —>■ 0. However, it is possible to combine 
the singularities in the forth line so that we find only one singular point and also 
an accompanying factor of two thus providing a cancellation of singularities in 
the third and forth lines. 

Indeed, the expression under the integral in the forth line in Eq.(|14D is sym- 
metric under the simultaneous transformation k\± k 2 ± and Ay «-> —Ay which 
is nothing else than the permutation of the two produced particles. Therefore we 
can simply multiply this term at 28(ki± > k 2 ±) to obtain: 

do r 



d 2 k±dy 

J d q 1± d q 2± - 2 - 2 J d k 2± dAy————— 

+2 J [ q 2 ± d 2 k ± d 2 k 2± dAy q 2 ± ^ **> > ^ ^ 



<p(xi,o,qi±) da 2 v(x 2 ,o, qa±) 



+ / d$ 



g 2 _i_ d 2 ki±d 2 k 2 ±dAy q 2j _ 
( <f(xi, <?i_l) </?(x 2 , g 2 ±) ^(^1,0, gi±) ^(^2,0, <?2± 



6{R{k 1 ,k 2 )>R )6{k 1± >k 2± ) 



Qi± qi± qf± Q2± 



f/ " 2 -0(J2< /,'„). 



d 2 ki±d 2 k 2 ±dAy 



The combination of the third and forth lines in Eq.(|T|) is free from singular- 
ities (note that for k 2 — ► there is no difference between R(k, k 2 ) and R(ki, k 2 )) 
and the singularity in the second line cancels with that in V in Eq.(^). Combining 
Eq.fllSf) with Eqs. (§),@ we obtain the second order correction to the one-particle 
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inclusive production in the high energy factorization scheme: 



d 2 k±dy 

{ d 2 q^d 2 q z ^ ^ Xlfi ' ^( X2 '°' q2± "> ( d ° 2 + d ° v \ 

o/^ M X l'^) ^2 <^(x 2 ,q 2± ) , , p v 



d 2 ki±d 2 k 2 ±dAy q% ± 

+ I d<S> 

-9(R < Rq). 



-6(R{k h k 2 )>Ro)9(k 1± >k 2± ) 



(p(xi,qi±) y(x 2 ,q 2± ) _ fjxifi, gi±) y(^2,o,92±) 



qI± ql± qf± q 2 2L 



d 2 k\±d 2 k 2 ±dAy 

Let us now analyse the MRK limit of Eq. (|i6|) , that is we take limit Ay — ► oo 
in ^4 (see Eq.(0) and Eq. (|A.2 ) in Appendix). In this limit 



A -> Amrk = X A (17) 

fi 'l± fi '2± 

which is precisely the combination of two leading order BFKL kernels responsible 
for real particles production. However, let us remind that the leading order one- 
particle production in high-energy factorization, described by Eq.(|) includes 
MRK contributions to all orders if the unintegrated structure function, <f>(x, q±) 
includes resummation to all orders of a s ln(l/x). It is evidently the case for 
structure functions undergoing BFKL equation. For other types of structure 
functions we just assume that the resummed a"ln n (l/x) terms are included in 
some hidden way. Consequently, we must subtract Amrk from A, and we will 
imply this subtraction in the following. 
To proceed further we must calculate 

d&2 - U^dAy- d&2 



d 2 k± J d 2 ki±d 2 k 2 ±dAy 

Note that the quantity in brackets in the first line of Eq. ([ITJ) should coincide 
(up to the constant factors depending on normalization) with the NLO BFKL 
kernel written explicitly in Note however, that when calculating the real 
contribution in JIB] the terms vanishing after integration over d 2 k± were dropped 
which did not change the result for NLO BFKL Pomeron itself. In calculating 
the one-particle inclusive cross sections these contributions have to be kept. 

The integration in Eq.flT8|) is very difficult. Fortunately, we can do it not for 
whole 02, but only for its singular part erf. We also have to change some other 
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terms in Eq.(|T6|) that emerge when arriving from Eq.(|i~2"|) at Eq. (|T6|) . Finally, the 
result is 



d 2 k±dy 

[ d 2 q L± _d 2 q^ V^ ^ 1 ' 0, gl± ) fl^Ml g2± ) ( d ° 2 + ^1 
J ah ah \d?k\ d 2 k 



Til T2± 

PtU^U/'T 1 ' ^'^' > *> (19) 

+ [d$( (p ( Xl,qi± ) ( P( X 2,Q2±) 

J \ q 2 ± d 2 k\±d 2 k 2 ±dAy q 2± 



Qi± d 2 k i± d 2 k 2± dAy q 2 ± 



6{R < Rq). 



3.2 Integration. 



Let us now choose the singular part of A with contributions from quark and gluon 
production added up in the form (see Appendix, MRK part as it was mentioned 
above is subtracted): 

A s = ^fi?L . gig2 2 ch(Ay) _ / _ nj \ 2g|g| f D-2 _nf_\^ , , 
2^2^ + fciA; 2 s v 4ivJ + I 2 NJ 8s 2 { ] 

This form of A s is chosen to avoid artificial ultraviolet divergency which occurs 
if one takes S = k 2 as it is in collinear and infrared limits. For the same purpose 
we take E in the form: 

E = | [^fe - (fci - fc) - (?? - g 2 2 ) (fc? - k 2 2 ) + 2hk 2 sh(Ay) (k 2 - q\ - q\)\ (21) 

Now we integrate the singular part of A over transverse two-dimensional mo- 
mentum space analytically continued to D — 2 = 2 + 2e or, strictly speaking, 
with 

„, d 2+2e k 2± , . 

and over rapidity. The results of the calculation of the integral over A s are given 
in Appendix B. The answer reads (see Eq.([|) for relation between a 2 and A) 

da s r 2N 2 a 2 6^(q 1± + g 2 ± - fcjj T(l - e) 4r 2 (l + e) ( k 2 ^ £ 

d 2 k ± ~ N 2 -l k 2 (47r) 1 + e er(l + 2e) \fi 2 

(- + 2^(1) -2^(1 + 2e) 



(23) 



1 . „ n ll + 8e n f A + Qe 



2(l + 2e)(3 + 2e) AN C (1 + 2e)(3 + 2e) 



The result for a v was derived in pOfl . Note that the answer depends on the 
arrangement of different corrections to QMRK amplitude. In this paper the 
symmetric variant ffl is chosen. 



da v 4N 2 a 2 s 6®{q 1± + g 2 ± - k±) T(l 




ql - ql 

<i. _ £t _ 2 in ^] 2q * q2 - MgM + q2 + A q^) + mA 
q\ ql n ql) (ql-ql) 3 qhl) 



(24) 



From Eqs. (|2"3"|) and (^) it is easy to see that divergencies of real and virtual 
parts cancel leaving a finite contribution to the first line in Eq . (|i~9l) : 




11 2nA, k 2 
In — 



3 3NJ (i 



2 



(25) 




Si _ M _ 2 In ^] 2 gigi ~ q\qi{ql + ql + 4gig 2 ) qiqA 

ql ql n ql) (ql-ql) 3 qhl). 

The first term in Eq. (p5|) , which is proportional to ln(/c 2 //z 2 ), is nothing but the 
well-known contribution corresponding to the running coupling. Therefore, after 
replacing a s by the running coupling a s (k 2 ), one should drop this term. 

The Eqs. ( |19|) and (p5|) together with Eq . (|20|) and the formula from Appendix 
A provide an analytical expression for the one-particle inclusive cross section. 

For practical applications one should integrate over the parameters of unin- 
tegrated structure functions (disregarding trivial elimination of delta-functions). 
The corresponding numerical calculations numerical studies will be descried in 
the next Section. 



4 Numerical estimates 



The numerical results strongly depend on the type of structure functions used 
in the calculation. Let us first consider the asymptotic BFKL structure function 



21 



(p(x,q 2 



x 



q 



qo ^Tc\"\n(x /x 



: exp 



ln 2 (g 2 /g, 



4A"ln(x /x) 



(26) 



with A = 4\n2N c a s /ii and A" = 14£ (3)N c a s /ii; a s is chosen to be equal 0.2. 
Because the asymptotic BFKL structure function is a solution of the linear ho- 
mogeneous equation, it does not have definite normalization and, moreover, the 
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parameters go and %o are arbitrary. Since the calculation with this structure func- 
tion is illustrative only, we simply choose C = 1, go — lGeV and Xq — 1. From 
Eq.(|T9|) it is clear that the NLO correction to the production process depends 
on the parameter R describing the collinear angle. For this calculation we take 
R = 0.7. Furhter discussion of the R dependence will be given below. 

It is important to note that although in asymptotic BFKL structure function 
the strong coupling constant does not run, to be consistent we should however 
make it run in a semihard vertices (cf. the discussion after Eq. (|T9|) ) . In the actual 
calculation with 1-loop a s we choose Aqcd = 200MeV. 

The one-particle inclusive cross section for v^S* = 5.5TeV, y — 0, and rif = 4 
is shown in Fig. 1 where for collinear factorization Eq.([I|) was used with 

dq 2 



xg(x,k 2 ) = J -^-(j)(x } q 2 ) 
o ^ 



As our second example we study the production process using the AKMS 
structure function |[22|| . In AKMS approach unintegrated structure function is 
obtained by application of BFKL equation to the evolution of structure function 
on x from x = 10~ 2 where <f)(lO~ 2 , k 2 ) is chosen to give the best fit for experi- 
mental data. 

This structure function can be fitted with an acceptable accuracy at least in 
the range of interest x = 10~ 2 ..10~ 4 by 



, 2\ A 9 
[x,q) = -g — exp 
xPq 



B(x)\n z (q 2 /q' 2 (x)) 
41n(l/x) 



(27) 



with parameters A = 8.55 ■ KT 2 , (3 = 0.486, q' = 0.758 + 5.41a; - 6816 and B = 
3.51 — 2.91 lnln(l/a;) + 0.793(lnln(l/s)) 2 . The results of calculation are given in 
Fig. 2 again with 1-loop a s , A QC d = 200MeV, VS = 5.5TeV, n f = 4, y = 0, and 
R = 0.7. 

Finally, in Fig. 3 we show the one-particle inclusive cross sections calcu- 
lated with GRV94(NLO) structure function P3| satisfying the DGLAP equation 



not related to BFKL. However, it possibly includes leading MRK contribution 
through the initial conditions for DGLAP evolution. Note the rapid increase and 
broadening of structure function with decreasing x, which is the characteristic 
property of BFKL induced structure functions. This calculation was done for 
= 1.8, 5.5 and 14TeV, 2-loop a s with A.q C d = 200MeV (because with these 
parameters GRV94(NLO) is calculated), rif — 4, y — 0, and R = 0.7. 

From Figs. 1,2,3 we see that taking into account the NLO corrections leads 
to the decrease of the particle production rate at high energies. For the structure 
functions obtained in BFKL approach NLO corrections change cross sections 
substantially (up to 50% in chosen kinematical interval), for non-BFKL GRV 
structure function changes are more dramatic: corrected cross sections are 2 to 
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5 times smaller than the leading order ones (apparently this results are sensitive 
to the cone size). We have no explanation to this fact. We hope that further 
studies of relationships between different types of structure functions may make 
this subject clearer. 

Figs. 3a,b,c show that ratios of differential cross sections calculated to the 
LO and NLO accuracy in high energy factorization and to the LO in collinear 
factorization are insensitive to the c.m.s. energy in the domain x < 1. 

Finally, in Fig. 4 we show the dependence of NLO cross section on the cone 
size R. It may be fitted well by function of the type A + B In R + CR and becomes 
infinitely large (and negative) at R — > 0. This is the general property of quantities 
with cancelling virtual and real corrections showing that at small values of R the 
fixed order perturbation theory is not valid (cf. ||24j| ). 

5 Conclusion 

The paper is devoted to the calculation of next-to-leading order correction to 
one particle inclusive production in the framework of high-energy factorization. 
High energy factorization scheme allows one to account for the initial transverse 
momentum of the colliding partons. The natural setup for particle production 
processes leading to high energy factorization is provided by Quasi-Multy Regge 
Kinematics. 

The results of computation of NLO contributions to BFKL Pomeron (cf . [[TJ 
and references therein) can be used to compute the next-to-leading order cor- 
rections to one particle inclusive production at high energies. This correction 
includes real and virtual pieces. The infrared singularity in the virtual piece in 
the NLO contribution cancels the infrared singularity in its real one when an 
infrared safe jet algorithm is applied. The explicit calculations of the infrared 
stable one-particle inclusive cross section at the next-to-leading order constitutes 
the main result of the paper. 

Numerical estimates were made to analyze the magnitude of NLO corrections 
for typical semihard transverse momenta and central rapidity region. Here we 
observe an essential dependence on the type of the structure functions used in 
the computation and shows more stability for BFKL-type structure functions and 
than for DGLAP one. 
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Appendix A Cross sections of pair production in 
high energy factorization 

We will use the following notation: 

s = 2(k 1 k 2 ch(Ay) - k 1± k 2± ); 

t = -(qi± ~ h±) 2 ~ hk 2 e Ay , u = -(q 1± - k 2A _) 2 - kik 2 e~ Ay ; 
£ = %\X 2 S = kf + k 2 + 2kik 2 ch.(Ay) , 



with k\ = Jki X , k 2 = yk 2X and ki±k 2 ± is the dot product with 2d Euclidean 
metric. 

Combined gluons and quarks (fermions) contribution to gg scattering has the 
form (adopted to fl5|)) 



<A. -A-gluons ~\~ ^jy3"^ fermions (^■•^-) 



c 



A.l gg ->- gg 



•Agluons *A\ ~\~ >A. 2 



22 f_l 1 q\q\ e Ay e~ Ay 1 
1_ qi<l2 1 tu + Atuk\k 2 2 Atk x k 2 Auk l k 2 + Ak\k 2 2 + 



1 

E 



- (l + k x k 2 {- ~ -)sh(Ay)) + -^(1 + -)ch(Ay)- 
s \ t u J lk\k 2 s 



4s ki t k 2 u 

As k 2 t kx u 



. D ~ 2 jf (ki ± - q 1± ) 2 (k 21 _ ~ qi±? - k 2 k 2 \ 2 

A *~ — {[ Tu J 

1 ( (k 2± -q 1± ) 2 -kik 2 e~ A y _ E\ ( (fcn - gn) 2 - k x k 2 e Ay E\ ' 
4 \(k 2± - qi ± ) 2 + hhe-^v s J I (fen. - qi ± ) 2 + hk 2 e^ s J 



(A.2) 



2 2 

£ = (<7i± - g 2 ±)(fen - h ± ) - ^(q't - ql)(k\ - k 2 2 ) + 2* 1 fesh(A 1 /) ( 1 - 91 
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A.2 gg -> 



A/ = {2 f( 1 + faM .(A»)(i - i)) - ( <^ ~ M t = 911)2 = ^ V + 



sE \ t u J \ tu 

1 / (fc 2 ±-gi±) 2 -fcifc 2 e- A « E\ ( (k 1± -q 1± ) 2 -k 1 k 2 e^ E\\ 

2 ^ (fc 2± - gi ± ) 2 + fcifee-^ s J ^ M - <M 2 + fcifc 2 e A » s)j 1 ' j 

and 

A/ = | | M±-^ 2 M-gixM^ 2 _ qjgj 
where E is the same as for gluons. 



Appendix B Integrals 

To make integration easily it is worth to change integration over Ay to integration 
over x = k\/{k\ + k 2 e Ay ). 



i 



— < B1 > 



x(l — x) 

v 



d 2+2e A; 2± MM 1 f d 2+2s k 2± 



(2tt) 2£ V 2 kfki J " 2 J (2tt) 2£ (k ± -k 2± ) 2 k 2 



2± 



T(l-e) k 2 ^ T 2 (l + e) 
(Any k 2 eF(l + 2e) 1 ' ' 



r d 2+2 "k 2± ch(Ay) r d 2+2 "k 2± (l-x x \ 

J (2vr) 2e hk 2 s J (2tt) 2£ \2k 2 2 sx 2k 2 s(l - x) ) [ ' j 

f d 2 ^k 2± (1-x) 2 

J (2n)* 2k 2 ((l~x)k ± -k 2± ) 2 + {X " X) 
T(l-e)k 2e r 2 (l + e) , 2£ 2el 

= "l^FM&j ,( ~ ^ ] 

Let's now combine these two contributions and perform an integration over x 

r(i- £ )** r 2 (i + E ) f dx „, , 2e , , e 



I » = '^T?7r^)J*r4 [{1 - x)+ ^- 1] - iBA) 
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In order to avoid divergencies we introduce infinitesimal parameter 8 (5 <C e) so 
that 



fdx 22£ _ = lim [ dx [(l-rr) 2e + x 2e - 1] 

and 

'» = "WTO<I + 2 * (1) " 2W + 2£)) (R6) 

For the integrations involving S it is useful to change k 2 ± on k — k 2 ± — (1 — 
x)k± so that s = k 2 /x(1 — x) and £ = /t 2 /x(l — x) + /c 2 . Now 

y dx r d 2+2e k 2± 1 = T(l-e) k 2£ T 2 (l+e) 1 
J x(l -x)J (2tt) 2£ sS ~ n (Any k 2 eT(l + 2e) 1 + 2e { ' ' 



and 



dx fd 2+2£ k 2 ± E 2 T[l-e)k 2£ T 2 (l + e) l-e 

i — 7]-. 



(2tt) 2£ 8g 2 g|s 2 (4vr) £ k 2 eF(l + 2e) 2(1 + 2e)(3 + 2s) 

(B.8) 
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Figure captions 

All cross sections da/cPkdy are calculated for combined gluon and quark con- 
tributions with n/ = 4, y — 

Figure 1 da/d 2 kdy calculated with asymptotic BFKL structure function Eq . (|26|) . 
VS = 5.5TeV, 1-loop a s with A QCD = 200MeV, R = 0.7 

Figure 2 da/d 2 kdy calculated with AKMS structure function Eq.(^). = 5.5TeV, 

1- loop a s with A Q cd = 200MeV, R = 0.7 

Figure 3 da jd 2 kdy calculated with GRV94(NLO) structure function [2l|. 2-loop a s 
with A QCD = 200MeV, R = 0.7 

a. = 14TeV 

b. VS = 5.5TeV 

c. = 1.8TeV 

Figure 4 .R-dependence of dojd 2 kdy calculated with GRV94(NLO). y/S = 5.5TeV, 

2- loop a s with A QCD = 200MeV. 
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Figure 4 
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